Quantum mechanics of spin transfer in coupled electron-spin chains 
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The manner in which spin-polarized electrons interact with a magnetized thin film is currently 
described by a semi-classical approach. This in turn provides our present understanding of the 
spin transfer, or spin torque phenomenon. However, spin is an intrinsically quantum mechanical 
quantity. Here, we make the first strides towards a fully quantum mechanical description of spin 
transfer through spin currents interacting with a Heisenberg-coupled spin chain. Because of quantum 
entanglement, this requires a formalism based on the density matrix approach. Our description 
illustrates how individual spins in the chain time-evolve as a result of spin transfer. 
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The problem of switching the spin state of a ferro- 
magnet has its origins in magnetization reversal through 
an applied magnetic field The present understand- 
ing of this process is based on the (numerical) solution 
of the so-called Landau-Lifshitz-Gilbert (LLG) equations 
0. These equations define a space-dependent magneti- 
zation vector, whose magnitude at each spatial coordi- 
nate is held fixed, but whose direction is determined by 
a (classical) dynamical equation which relates the rate of 
change of magnetization to the applied torque. The ap- 
plied torque is determined by an effective magnetic field, 
which represents a mean field as experienced by the mag- 
netization at that particular spatial coordinate. A damp- 
ing term is also required, which is crucial for the magne- 
tization reversal to take place. A semi-classical approach 
[H, Q was adopted to describe spin transfer [j| from a 
spin-polarized current to a ferromagnet. The descrip- 
tion of this phenomenon requires a quantum mechanical 
treatment of the spin current, while a classical notion of 
the magnetic moment [H 0, Q suffices for a qualitative 
understanding of the spin torque process. Most intrigu- 
ing from our point of view is that a phenomenological 
damping process is not required for spin torque to occur; 
in this respect a spin current is profoundly different from 
an applied magnetic field. 

However, more recently, spin transfer has been recog- 
nized as important, both in spintronics 
[HI], and in quantum information (Til, fljj]. For these rea- 
sons, and because, after all, spin is an intrinsically quan- 
tum property, it is critical to understand spin transfer 
as a fully quantum mechanical process. This way, even 
within the more conventional applications where the LLG 
approach has been fairly successful, one can understand 
and perhaps push the limitations of the semi-classical 
approach. In this paper we construct the simplest but 
most important elements of a model with which one can 
illustrate the spin dynamics quantum mechanically. Con- 
sider an electron incident on a ferromagnetic chain with 
N s coupled local spins, all prepared magnetized in a par- 



ticular direction. Unless the electron spin is parallel to 
the local spins, spin will be transferred to the chain from 
the electron. In the semi-classical picture, each spin (S) 
is considered as a classical vector; in particular the mag- 
nitude of the vector remains constant, while its direction 
is allowed to change. On the other hand, within a quan- 
tum mechanical description, S is an operator USES, it 
is (S 2 ) that remains constant; in contrast, (S) changes 
during the transfer process, and can even vanish momen- 
tarily 

To describe complete spin transfer requires a current 
of polarized electrons. This introduces another purely 
quantum mechanical concept associated with spin trans- 
fer: entanglement between an electron and the localized 
spin chain. Even after the electron moves away from the 
spin chain, the quantum states of the electron and local- 
ized spins are entangled in spin space. Thus, to illustrate 
the dynamics of the chain we should know how to de- 
scribe a subsystem (spin chain) of a composite system 
(spin chain and electron). In particular, if we send an- 
other electron to interact with the chain after a preceding 
electron leaves the chain, it is crucial to decouple the spin 
chain from the first electron quantum mechanically. The 
density matrix formalism [1 81 ] is the only way to serve 
this purpose. Suppose a system consists of two subsys- 
tems A and B. This composite system is assumed to be 
closed. Initially, A and B are separated and no interac- 
tion takes place between the two. In other words, A and 
B are further apart than the interaction range of the two. 
The quantum states of A and B would be \ipji) and \i/jb), 
respectively. Imagine, now, the two subsystems come to- 
gether and interact with one another. For example, in 
our case, the mobile subsystem A (the electron) moves 
into the region of interaction with the stationary subsys- 
tem B (the localized spin chain), they interact, and then 
A leaves the domain of B . Let us assume the interaction 
lasts for a finite amount of time, say tQ. The dynamics 
of the composite system will be described by the total 
Hamiltonian H = Ha + Hb + Hi nt , where Ha and Hb 
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are the Hamiltonians of A and B, respectively, while Hi nt 
is the interaction between the two. 

The density matrix p tot of the total system at t = 
is p tot (0) = p A (0) ® p s (0), where p A (0) = |^>(^a| and 
p B (0) = \i^b)(iPb\- The quantum states of A and B 
are assumed to be pure initially, which means that these 
states are exactly known at t — 0. The time evolution of 

pt0t j g gj V6n 



P tot (t) = U(t)p tot {0)uHt) 



(1) 



where U(t) — e lHt . The expectation value of an opera- 
tor is calculated as follows: 

(SI) = T^a,b W ot {t)} = Y,(*\mp tot {t)\b)\a) (2) 

a , b 

where {|a}} and {\b)} are orthonormal sets associated 
with the subsystems A and B, respectively. In partic- 
ular, for an operator SIb of B, the expectation value is 
(SIb) = Tr# [£!bp b ] where the reduced density matrix 
of the subsystem B, p B , is obtained by tracing out the 
states of the subsystem A: p B = Tr A [p tot (t)]. While this 
relation is obvious, it serves to indicate the significance of 
the reduced density matrix p B . After the two subsystems 
no longer interact at t = to, the quantum states of B are 
represented solely by p B . The quantum state of B is not 
pure but mixed; in other words, it is impossible to repre- 
sent p B (to) as |4 r s)(^ r s|, where \^b) is a quantum state 
exactly known at to- It is a generic property of quantum 
entanglement that a subsystem is in mixed states after 
getting decoupled even if the composite state is pure [14 1. 
If another subsystem (second electron) is introduced and 
interacts with B at time t , p B (t ) should be used as a 
part of the initial density matrix corresponding to sub- 
system B. Subsystem A can be introduced as for the first 
electron, but representing the next electron. Then the 
initial density matrix is p tot (to) — p A (to) ® p B (to), where 
P A (to) — P A {0) by preparation, and its time evolution 
will be p tot {t) = U(t) [p A (t Q )® p B (t )] W{t) for t > t . 
We then apply the same procedure systematically in this 
way to allow a beam of electrons to individually interact 
with the spin chain. 

As a model Hamiltonian [h| for spin transfer from 
incident electrons to a localized spin chain we use 
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1 = 1 



1=1 



(3) 

where c ia creates an electron with spin a at site i, S/ is 
a localized spin operator at site I, t is a hopping ampli- 
tude between nearest neighbor sites, so that the electron 
can move along the length of the entire chain, and Jo is 
the coupling between the electron spin and a local spin. 
Note that Jo coupling takes place only when an electron 
is on the site of the local spin. The parameter J\ is a 
coupling between two neighboring localized spins. For a 



ferromagnetic chain J\ > 0. However, this formalism can 
be applied for arbitrary values of Jo and Ji; for example, 
one can study spin transfer to an antiferromagnet [201 ] 
using this model Hamiltonian. Moreover, these coupling 
parameters can be site-dependent. A convenient notation 
is to represent Jo as a vector with N s components and Ji 
as a vector with N s — 1 components whenever necessary. 
Since we do not include longer range spin-spin interac- 
tions, which would lead to domain formation, we can use 
a site-dependent J\ to mimic domains. Scalar potentials 
are not taken into account because they are irrelevant to 
spin transfer. 

Our calculation scheme is as follows: we propagate 
electrons, one after another, towards the ferromagnetic 
spin chain. These electrons are spin-polarized in the +z 
direction, while, initially, all local spins are aligned (for 
example, along the x, or the — z direction). Depending 
on the momentum of the incident electron, a time to can 
be determined, after which the electron and the spins no 
longer interact. At to, the reduced density matrix of the 
spin chain is evaluated from the total density matrix by 
decoupling the spins from the electron. At the same time 
we introduce another electron to interact with the spins. 
This reduced spin density matrix is used to construct an 
initial density matrix of the total system (spin chain and 
new electron). This process is repeated until eventually 
all local spins become aligned with the incident electron 
spin. 

The density matrix of the total system (spin chain 
and electron) is p tot (0) = p s (0) ® p el (0) at t = 0, 
where p el {0) — |0o)(<fo|- The initial state of the elec- 
tron is defined as |0o) = '^2 i tPi\h+) , where the am- 
plitudes ipi are chosen to represent a Gaussian wave 
packet with a mean momentum ko and a mean posi- 
tion Xq. We choose fc = 7r/2, and x sufficiently far 
from the chain such that initially the wave packet is well 
defined outside the chain. The spin state of the inci- 
dent electron is |+), i.e. along the +z direction. The 
reduced density matrix of the spin chain can be obtained 
by quantum-mechanically decoupling the chain from the 
electron: p s (t) = J2i o-ih °1/°' 'WK' cr )' where \i,a) is 
a state vector of an electron with spin a at a site i. 
As mentioned, the time evolution of p tot is given by 
p tot (t) = Up to \0)U\ where U = exp[~iHt}. At t we in- 
troduce another electron to interact with the chain. Now 
the initial density matrix of the spin chain and this elec- 
tron is p tot (t ) = p S {to) ® p e/ (0), and its time evolution 
is p tot (t > t ) = Up s {t )®p el (0)W. After the (A + l)th 
electron, we obtain 



P S N+l(t) 



(i N ,<rN\Upff (Nt )U^\i N ,a N ) . (4) 



In our numerical calculation, we consider a ferromagnetic 
chain in a state \F) with N s local spins; thus, p (0) = 
\F)(F\. However, any initial spin configuration can be 
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FIG. 1: Time evolution of 6 spins interacting with 4 in- 
coming electrons one by one. The site-independent electron- 
spin coupling is Jo = 2 and the ferromagnetic coupling is 
Ji = 0.4, which is also site-independent. The arrows depict 
(S) as three-dimensional vectors for spatial as well as tempo- 
ral behavior of individual spins, and (S z ) is further visualized 
by colors. 



assumed. The expectation value of a spin operator, say 
S u , is evaluated by (S lz ) = Tr s [S lz p s ] . 

In order to describe spin transfer from spin-polarized 
incoming electrons to a spin chain, we show dynamics of 
the spins in the chain. We consider 6 spins with S =1/2 
initially pointing to the x direction. For the parameters 
we choose, the electron wave packet no longer interacts 
with the chain at to = 60 (in dimensionless units). Fig. 1 
illustrates how the 6 spins time-evolve while interacting 
with 4 incoming electrons one after another. The site- 
independent electron-spin coupling is Jo = 2 and the 
ferromagnetic coupling is J\ — 0.4, which is also site- 
independent, in units of the hopping parameter t(= 1). 
The arrows depict (S) as three-dimensional vectors for 
spatial as well as temporal behavior of individual spins, 
and (S z ) is further visualized by colors. As discussed in 
Ref.[17j|, the coupling J determines the degree of spin 
transfer, which depends on Jo non-monotonically. Fig. 1 
shows two distinct dynamical regimes. First, while an 
electron wave packet is striking the spin chain, the Jo 
interaction term in H causes spin transfer, and therefore 
spin flip. This will occur with every member of the spins 
with multiple reflections of the electron within the chain 
included. Once the electron has left the spin chain, how- 
ever, the spins continue to evolve in time, but solely due 
to the Ji term. The result is the oscillating behavior 
depending on the value of Ji . 



FIG. 2: Time evolution of 6 spins with a domain structure 
formed by Ji interacting with 4 incoming electrons one by 
one. The coupling Jo remains unchanged as in Fig. 1 while 
the site-dependent Ji = (5, 0.4, 5, 0.4, 5) for 3 domains. 



As mentioned earlier, parameterizing components of 
Ji, one can envision a domain-structured spin chain. In 
Fig. 2, the coupling Jo remains unchanged as in Fig. 1 
while we choose a site-dependent Ji = (5, 0.4, 5, 0.4, 5) 
for 3 domains to form. The dynamics of spins is very 
different from the previous case. Because of the strong 
coupling between three spin pairs, two spins in each pair 
evolve simultaneously. The first and the last spin pairs 
oscillate against each other while the middle pair evolves 
almost monotonically. This overall behavior of three spin 
pairs is qualitatively the same as one of weakly coupled 
three local spins. Note that as more electrons are sent 
very little changes; for these parameters, five electrons 
are nearly sufficient to align the entire chain along the 
+z direction. 

Another issue we wish to address in this paper con- 
cerns the excitations in the spin chain induced by the 
spin transfer process. We do not include a magnetic field 
in the Hamiltonian Eq. ^ because we investigate exci- 
tations only due to the relative spin configuration in the 
spin chain. If a magnetic field is included, additional exci- 
tations will be induced depending on the orientation of lo- 
cal spins with respect to the field. For investigation of the 
induced excitations, we consider only site-independent 
couplings Jo and J\. The excitations induced by spin 
transfer are, in fact, associated with the energy transfer 
between the incident electrons and the spin chain. The 
energy of the chain is given by the expectation value Es 
of the Hamiltonian of the chain H$ = — 2Ji J^i ■ Sj+i; 
namely, Es — Trs [Hsp s ] ■ We demonstrate that Es 
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FIG. 3: Excitation induced by spin transfer Es normalized 
by Es(0) = Ji (N s — 1) /2 as a function of the number of 
electrons sent in. The spin chain interacts with 25 incident 
electrons with Jo = 1 and various values of J\. The momen- 
tum ka of incident electrons is 7r/2 in the main frame while it 
varies in the inset. 



of spin transfer; spins as quantum operators and entan- 
glement. 

As far as the physics of spin transfer is concerned, 
it is clear that the semi-classical assumption that the 
magnitude of the magnetization vector remains constant 
during spin transfer is generally violated. In particular, 
when the ferromagnetic exchange coupling is weak we 
find regimes where the magnitude of the spin vector can 
almost vanish. How this might alter our understanding of 
experimental observations will be the subject of further 
investigation. 

This work was supported in part by the Natural 
Sciences and Engineering Research Council of Canada 
(NSERC), by ICORE (Alberta), and by the Canadian 
Institute for Advanced Research (CIAR). Numerical cal- 
culations have been performed using computational fa- 
cilities of the Academic Information and Communication 
Technologies (AICT) and the WestGrid. 



depends on the coupling parameter J\, and energy is 
not monotonically transferred to the chain from incoming 
electrons. As local spins start to align with one another in 
the direction of the incident electron spin, energy would 
be transferred back to incoming electrons from the chain 
as illustrated in Fig. 3. This backwards energy trans- 
fer can be understood in a sense that excitation and de- 
excitation are relative in the absence of a magnetic field. 
Fig. 3 is for a spin chain interacting with 25 incident elec- 
trons one by one. We set Jq = 1 but choose various values 
of J\. The mean momentum fco of incident electrons is 
7r/2 in the main frame while it varies in the inset to show 
that the energy transfer also depends on the momen- 
tum of the incident electrons. If J\ is large enough, say 
Ji > 5, then no excitations occur and Es/Es(0) = — 1, 
where Es(0) = J\ (N s — 1) /2 for a chain with Ns local 
spins. Since there is no preferred direction due to the 
absence of a magnetic field, Es remains unchanged as all 
spins rotate in unison. This can be shown analytically. 

In summary, we have outlined a quantum-mechanically 
rigorous framework with which one can calculate the spin 
transfer from a spin-polarized electron current to a spin 
chain. A brute force method is possible, but for a pro- 
gression of incident electrons, this will quickly expand 
the Hilbert space to exceed the reach of modern comput- 
ers. Application of the density matrix formalism not only 
solves this difficulty but also correctly deals with quan- 
tum entanglement inherent to the spin transfer problem. 
We have performed computations for a rather simplistic 
system; nonetheless, included are two essential elements 
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